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O I Abstract 

Q_i' In this paper the formulae are cohected which are needed for the computation of 

^ . the strong couphng constant and quark masses at different energy scales and for 

QJ [ different number of active flavours. All equations contain the state-of-the-art QCD 

' corrections up to three- and sometimes even four-loop order. For the practical 
implementation Mathematica is used and a package containing useful procedures is 

^ ! provided. 
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Program Summary 

Title of program: RunDec 
Available from: 

Ettp : //www-ttp . physik . uni-karlsruhe . de/Progdata/ttpOO/ttpOO-05/ 

Computer for which the program is designed and others on which it is operable: Any 
work-station or PC where Mathematica is running. 

Operating system or monitor under which the program has been tested: UNIX, Math- 
ematica 4.0 

No. of bytes in distributed program including test data etc.: 65000 
Distribution format: ASCII 

Keywords: Quantum Chromodynamics, running coupling constant, running quark 
mass, on-shell mass, MS mass, decoupling of heavy particles 

Nature of physical problem: The values for the coupling constant of Quantum Chro- 
modynamics, actually depends on the considered energy scale, fx, and 
the number of active quark flavours, rif. The same applies to light quark masses, 

m^^^\n), if they are, e.g., evaluated in the MS scheme. In the program RunDec all 
relevant formulae are collected and various procedures are provided which allow for 
a convenient evaluation oia^J^-^\fj) and niq'^\jj) using the state-of-the-art correction 
terms. 

Method of solution: RunDec uses Mathematica functions to perform the different 
mathematical operations. 

Restrictions on the complexity of the problem: It could be that for an unphysical 
choice of the input parameters the results are nonsensical. 

Typical running time: For all operations the running time does not exceed a few 
seconds. 

1 Introduction 

Quantum Chromodynamics (QCD) is nowadays well established as the theory of strong 
interaction within the Standard Model of elementary particle physics. In recent years 
there has been a wealth of theoretical results (for a review see At the same time 
perturbative QCD has been extremely successful in describing the experimental data with 
high precision. 

The fundamental quantity of QCD is the so-called beta function which connects the 
value of the strong coupling constant, at different energy scales /i. It is thus 
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particularly important to know the beta function as precise as possible. In the four-loop 
corrections were evaluated allowing for a consistent running at order a^. In the majority 
of all computations performed in QCD the MS renormalization scheme |^ is adopted. In 
this scheme the Appelquist-Carazzone decoupling theorem Q is not directly applicable. 
When crossing flavour thresholds, it is thus important to perform the decoupling "by 
hand". In order to be consistent, four-loop running must go along with the three-loop 
decoupling relation which was evaluated in p. 

Similar considerations are also valid for quark masses. Also here the renormalization 
group function is available up to the four-loop level [|] and the corresponding decoupling 
relation up to order ^ (see also P). 

In this paper all relevant formulae are collected which are necessary for the running 
and decoupling of as and for quark masses. Their proper use is discussed and easy-to-use 
Mathematica procedures collected in the package RunDec are provided. Their handling 
is described and examples are given. 

The outline of the paper is as follows. In the next Section the formulae are presented 
which are needed for the running of the strong coupling constant up to the four-loop 
level. The corresponding equations for the quark masses are presented in Section ^. In 
addition the conversion formulae between the MS and on-shell scheme are discussed in 
some detail. Section ^ is concerned with the decoupling of the strong coupling and quark 
masses. Finally, in Section |^, the most important procedures of the package RunDec are 
described in an easy-to-use way. For most practical applications they should be sufficient. 
In the Appendix the complete collection of procedures is given. 



2 Running strong coupling constant 

The beta function governing the running of the coupling constant of QCD is defined 
through 
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where Uf is the number of active fiavours. The coefficients are given by [Tl|, [T^, 
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( is Riemann's zeta function, with values (2 = and ^3 ~ 1.202 057. It is convenient 
to introduce the following notation: 
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In the following the labels fi and Uf are omitted if confusion is impossible. 
Integrating Eq. (|l]) leads to 
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where an expansion in a has been performed. The integration constant is conveniently split 
into A, the so-called asymptotic scale parameter, and C. The conventional MS definition of 
A, which we shall adopt in the following, corresponds to choosing C = (bi/Po) In Po R . 
Iteratively solving Eq. (Q) yields p 
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where L = ln(yU^/A^) and terms of 0{1/L^) have been neglected. A is defined in such a 
way that Eq. (|D does not contain a term proportional to (const. /L^) 0. 

The canonical way to compute a{^2) when a(/ii) is given for a fixed number of flavours 
is as follows: 



1. Determine A. There are several possibilities to do this. One could, e.g., use the 
explicit solution given in Eq. (^). Another possibility is the use of © and solve the 
equation iteratively for A. Furthermore the first line of (^) could be used and the 
integral could be solved numerically without performing any expansion in as- We 
will see in the examples below that the numerical differences are small. 

2. a(/i2) is computed with the help of Eq. where the value of A is inserted and /x is 
set to /i2. 

It is also possible to avoid the introduction of A in intermediate steps and to solve the 
differential equation (|I]) numerically using a(/i)|^=^^ = a(/ii) as initial condition. This 
convention requires the knowledge of both as and the scale /i in order to determine as at 
the new scale. Frequently, = Mz is used as reference scale. 
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On the other hand A plays the role of an universal parameter which at the same time 
sets the characteristic scale of QCD. 

At this point it is instructive to consider an example. Let us assume that as is given 
at the Z-boson scale: af\Mz) = 0.118. Let us further assume that it is determined from 
the experiment with three-loop accuracy, which means that in the (3 function (|^) only the 
coefficients up to [32 are considered and {3^ is neglected. Let us now evaluate the strong 
coupling at the scale /i = Mf, and compare the results obtained with the different strategies 
outlined above. In the following a possible Mathematica session is shown. NumDef is a set 
of Mathematica rules which assigns typical values to the physical parameters used in our 
procedures. The numbers used in this paper can be found in Eq. ([36| ) and the procedures 
are described in the Appendix. 

In[l] := «RunDec.ni; 

Comment: evaluation of A from af\Mz) based on the explicit solution, Eq. (^, and 
subsequent evaluation of af^M^) from A based on Eq. (^). 

In [2]:= lamex = LamExpl [asMz/ . NumDef ,Mz/ .NumDef , 5 , 3] 
Out [2]= 0.208905 

In [3] : = AlphasLam [lamex , Mb/ . NumDef ,5,3] 
Out [3]= 0.216610 

Comment: evaluation of A from af\Mz) based on Eq. (H), and subsequent evaluation of 
af\Mi,) from A based on Eq. (|]). 

In [4]:= lamim = Lamlmpl [asMz/ . NumDef , Mz/ . NumDef , 5 , 3] 
Out [4]= 0.208348 

In [5] : = AlphasLam [lamim , Mb/ . NumDef ,5,3] 
Out [5]= 0.216444 

Comment: evaluation of a^g\Mi,) from af\Mz) based on Eq. 

In [6] : = AlphasExact [asMz/ . NumDef , Mz/ . NumDef , Mb/ . NumDef ,5,3] 

Out [6]= 0.216712 

Rounding to three significant digits leads to a difference of ±1 in the last digit. Considering 
the direct integration of (|I]) as the most precise one we can conclude af\Mb) = 0.217 
assuming three-loop accuracy. 
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0.2059 


0.2059 


0.2059 


0.3369 


2 
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0.2173 


0.2161 


0.3965 


3 
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0.2164 


0.2167 


0.4029 
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0.2174 


0.2173 


0.2169 


0.4065 



Table 1: af\Mb) and ^(^^(l GeV) computed from af\Mz) usiner different number of 
loops for the running. For the computation LamExpl [] (a), Lamlmpl [] (b) (both in 
combination with AlphasLam [] ) and AlphasExact [] (c) is used. 



In Tab. |I| the influence of the number of loops is studied in the evaluation of af* at 
the scale and the (hypothetical) scale 1 GeV using af\Mz) as input. (For the latter 
only the function AlphasExact [] is used for the computation.) It can be seen that the 
inclusion of /3i leads to a significant jump in af\Mi,) whereas the effect of the three- and 
four-loop coefficients, i.e. j32 and /^s, is only marginal. Their influence is more pronounced 
for /i = 1 GeV. 



3 Quark masses in the MS and on-shell scheme 

In the MS scheme the running of the quark masses is governed by the function 7m(c^s) 
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where the coefficients 7m,i are known up to the four-loop order |Tj, |T5|, |T6 
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with Ca ~ 1.202 057, (4 = 7r790 and Cs ~ 1.036 928. In analogy to (D we define 



Combining Eqs. (|T]) and (H) leads to a differential equation for m(/x) as a function of 
asifi). It has the solution |T7 



with 
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where terms of 0(x'^) have been neglected. For a given mass, m, at scale /iq and a<((/Uo) 
the scale invariant mass /i^ = rn{fim) can be obtained from Eq. (P) by iteration. Note the 
appearance of as{fi) on the r.h.s. of (^. Thus for the computation of fim it is convenient to 
use in a first step Q;s(yUo) in combination with Eq. (H) to determine A. Afterwards Eq. (P) 
is used again for the calculation of as(/i) which is inserted in (||) before the iteration. 
From Eq. (P) it appears natural to define the mass 0] 
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which is often used in the context of lattice calculations. By construction the mass m 
is scale independent. It is furthermore scheme independent (as far as mass-independent 
schemes are concerned). This can be seen by considering the r.h.s. of (^^) in the limit 
II —>■ 00 
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and by recalling the fact that the coefficients /5o and 7^,0 are scheme independent. In the 
following we will refer to m as renormalization group invariant mass. 

In the following we want to provide the relations between the MS and the on-shell mass. 
Whereas the coefficient of order has been available since quite some time |l^, |19| only 
recently the three-loop result could be obtained |^D|, In an asymptotic expansion 
in combination with conformal mapping and Fade approximation has been used in order to 
obtain a numerical result for the MS-on-shell conversion formula. The numerical results 
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of 1^ are in perfect agreement with the subsequent analytical calculation of (cf. 
Tab. H). For a given on-shell mass the MS quantity can be computed with the help of 
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where (2 = 7r^/6 and I^m = In/x^/M^. is the number of light quarks. The function 
A(x) arises from the two-loop diagram with a second fermion-loop [l^. For < x < 1 it 
is approximated within an accuracy of 1% by 
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The corresponding mass effects at order are not yet known. In the argument of A(x) 
the ratio of the on-shell mass of the light quarks, Mj, and the heavy one, M, appears. 

The coefficients z^^\M) can be found in Tab. § for different values of ni where both 
the results of and pT| are listed. For completeness also the corresponding two-loop 
coefficients (without the contribution from A(x)) are given. The analytical result for 
zl^\M) reads [BT] 
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where 04 = Li4(l/2) ^ 0.517479. 

Iterating (^) leads to a relation between the scale-invariant mass, fim = rn{fim), and 
the on-shell mass 
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Table 2: Two- and three- loop coefficients of the relation between on-shell and MS mass. 
The choice fi"^ = M^, respectively, /i^ = has been adopted. 
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Inverting Eq. (|T3|) leads to 
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where for convenience /i^ = has been chosen. The numerical values of the coefficients 
and z^^'" can also be found in Tab. ^. Their analytic expressions are easily obtained 
from Eqs. (|13D and (0). 

Eq. ( p!7D can be used to compute the on-shell quark mass if the corresponding mass 
in the MS scheme is provided. In order to avoid large logarithms it is suggestive to use 
in a first step the renormalization group equation (^ and evaluate fim- In a second step 
Eq. (|r^) is used for fi = fim- Also in the case when the on-shell mass is given it is 
advantageous to use Eq. (|r^) for the computation of the MS mass. The reason is that 
Eq. ([T3|) contains contributions from the ill-defined pole mass of the light quarks like, e.g., 
the strange quark. In the case of the top quark it is safe to use (|T^) as in general the 
contributions for the charm and strange quark masses can be neglected. 

Concerning the determination of the quark masses a crucial role is played by lattice 
calculations. There it is not possible to use directly the MS scheme as it is tightly 
connected to dimensional regularization. Rather one has to use a prescription which is 
based on the so-called momentum subtraction scheme. In general these schemes have the 
disadvantages that they are not mass independent. Recently, however, a mass definition 
based on momentum subtraction — the regularization invariant (RI) mass — has been 



proposed which enjoys this feature In the relation to the MS mass has been 

evaluated to three-loop accuracy. It reads: 
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Table 3: Computation of ml '{Mz), Mj, and m^\Mz) from /if, = ml (/if,) = 3.97 GeV 
for different number of loops. 
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Let us at this point consider an explicit example. For a given mass 
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3.97 GeV five-flavour running (cf. Eqs. (^ and (p!0[) ) is used in order to obtain m\ \Mz) 
and the on-shell mass is computed with the help of Eq. (p!?]). Furthermore the value for 



m 



RI 



(Mz) is evaluated. In Tab. ^ the results are listed for different number of loops. 



At the end of this section we want to summarize the different mass definitions intro- 
duced in this section in the following table: 




on-shell mass 
MS mass 

scale invariant mass 
renormalization group invariant mass 
regularization invariant mass 



4 Decoupling at flavour thresholds 

In MS-like renormalization schemes, the Appelquist-Carazzone decoupling theorem Q| 
does not in general apply to quantities that do not represent physical observables, such as 
beta functions or coupling constants, i.e., quarks with masses much larger than the consid- 
ered energy scale do not automatically decouple. The standard procedure to circumvent 
this problem is to render decoupling explicit by using the language of effective field theory. 
The formulae presented below are valid for QCD with rii = nj — 1 massless quark flavours 
and one heavy flavour h, with mass m^ which is supposed to be much larger than the 
energy scale. Then, one constructs an effective n^-flavour theory by requiring consistency 
with the full n^-flavour theory at an energy scale comparable to m/j, the heavy-quark 
threshold /x^"-'') = 0{mh). This leads to a nontrivial matching condition between the 
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couplings and light masses, rrig, of the two theories. Although, a 
and m^^^\mh) = rri'q^\mh) at leading and next-to-leading order, this relation does not 
generally hold at higher orders in the MS scheme. At 0{a^g) the corresponding correction 
terms have been computed in ^ 

The connection between the strong coupling constant in the effective and the full 
theory is given by 



(19) 



where Q is known up to the three-loop order P, |^ : 
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In this equation the MS mass mh{lj) — indicated by the superscript MS — is chosen 
for the parameterization of the heavy quark mass and represents the renormalization 
scale. Often it is convenient to express C,g through the scale invariant mass, denoted by 
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Transforming the heavy quark mass into the on-shell scheme leads to 
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In practical applications also the inverted formulae are needed which read for 
Eqs. dO), dll) and (H): 
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The decoupling relations (^U|)-(^) have to be applied whenever a flavour threshold is to 
be crossed. 

At this point we briefly want to comment on the order of which has to be used 
for the running, respectively, the decoupling if the analysis should be consistent. If the fj, 
evolution of a^J^^\fi) is to be performed at + 1 loops, i.e., with the highest coefficient in 
Eq. ([1|) being [3^^'\ then consistency requires the matching conditions to be implemented 
in terms of A^-loop formulae. Then, the residual /i dependence of physical observables 
will be of order + 2. 

As an example let us compute a^g\Mc) from af\Mz) = 0.118. Let us furthermore 
consider the on-shell definition of the heavy quark, M^, i.e. we use Eq. (^) for the analysis. 
For the scale // in ([19D where the matching is performed we choose ^th = Mh. Assuming 
four-loop accuracy for the beta function and (as a consequence) three-loop accuracy for 
the matching computation one would proceed as follows (see Appendix for a description 
of the procedures): 



In [2]:= (alsmuth = AlphasExact [asMz/ . NumDef , Mz/ . NumDef , Mb/ . NumDef , 5 , 4] ) 
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number of loops 








(running) 


{fith = Mz) 


{f^th = Mb) 


ifith = 1 GeV) 


1 


0.3213 


0.2918 


0.2784 


2 


0.3387 


0.3318 


0.3222 


3 


0.3396 


0.3364 


0.3319 


4 


0.3399 


0.3378 


0.3350 



Table 4: Computation of a^J^\Mc) from a^^\Mz) = 0.118 for different number of loops. 
For the matching scale Mz, Mb and 1 GeV has been chosen. 



Out [2]= 0.2169467 

In[3]:= (alsmuthp = DecAsDownOS[alsmuth,Mb/.NumDef ,Mb/.NumDef ,4,4] ) 
Out [3]= 0.2163396 

In [4]:= (alsMc = AlphasExact [alsmuthp,Mb/ . NumDef ,Mc/ .NumDef ,4,4] ) 
Out [4]= 0.337848 

Finally one arrives at ai^^Mj = 0.338. These steps are summarized in the function 
AsRunDec [] where the corresponding call would read 

In [5] : = AsRunDec [asMz/ . NumDef , Mz/ . NumDef , Mc/ . NumDef , 4] 
Out [5]= 0.337848 

Note that the loop-argument of the function DecAsDownOS [] (last argument) refers to the 
order used for the running, i.e. in the considered case the "4" means that the three-loop 
relation is used for the decoupling. In this example the effect of the decoupling is quite 
small. It is actually comparable to the uncertainty from using different methods for the 
running (cf. Tab. |I|). However, one has to remember that for the matching scale the heavy 
quark mass itself has been used, whence all logarithms in Eq. (p^) vanish. A different 
choice would lead to a different result for a[^\Mc). On the other hand, on general grounds, 
the decoupling procedure should not depend on the choice of that scale, respectively, the 
dependence should become weaker when going to higher orders. In Tab. § the dependence 
of a^'^^Mc) on the number of loops is shown. For the matching scale Mz, Mb and 1 GeV 
has been chosen. It can be clearly seen that the four-loop analysis provides the most 
stable values — even in the case when the matching is performed at the a high scale like 
the Z boson mass. This is expected on general grounds as physical results should not 
depend on the matching scale. 

We should mention that in case the MS definition for the heavy quark is used in a first 
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step mh{fith) has to be evaluated. The corresponding formulae can be found in Section]^. 
They are also implemented as Mathematica procedures and described in the Appendix. 

In Fig. |l| it is demonstrated that the inclusion of the four-loop coefficient jSs accompa- 
nied by the three-loop matching leads to an independence of ^th = yU*'^'' over a very broad 
range [0. The plot shows the dependence of af\Mz) on the matching scale (denoted by 
^(5)) where a^'^^M^-) is used as starting point. Our procedure to get the different curves 
is as follows. We first calculate by exactly integrating Eq. (|lD with the initial 

condition a^f^M^) = 0.36, then obtain a^^H/^^^^) from Eqs. (gSj) with = 4.7 GeV, and 
finally compute a^^\Mz) with Eq. ([l|). For consistency, A^-loop evolution must be accom- 
panied by {N — l)-loop matching, i.e. if we omit terms of 0{a^~^^) on the right-hand side 
of Eq. ([5), we need to discard those of 0{a^) in Eq. (|25|) at the same time. In Fig. ||, the 
variation of a^^\Mz) with n^^^ /Mf, is displayed for the various levels of accuracy, ranging 
from one-loop to four-loop evolution. For illustration, fx^^^ is varied by almost two orders 
of magnitude. While the leading-order result exhibits a strong logarithmic behaviour, it 
stabilizes as we go to higher orders. The four-loop curve is almost fiat for n^^^ ^ 1 GeV. 
Besides the /x*-^-' dependence of af'\Mz), also its absolute normahzation is significantly 
affected by the higher orders. At the central matching scale /x^^-* = M^, we encounter a 
rapid, monotonic convergence behaviour. 

Fig. |l| can immediately be reproduced with the help of the procedure A1L2A1H[] 
described in the Appendix. 

Up to now only the decoupling of the coupling constant has been considered. However, 
also the (relatively) lighter quark masses undergo a decoupling procedure when crossing 
a fiavour threshold. If we define the connection between the quark mass in the effective 
and full theory through 



m 



(nf-l) 



(26) 



the decoupling constant Cm is given by 
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Note that all three quantities in Eq. 



depend on the renormalization scale /i. 
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(5) 




(5) 



Figure 1: fi^^'^ dependence of a'f\Mz) calculated from af\Mr) = 0.36 and Mb = 4.7 GeV 
using Eq. (|l|) at one (dotted), two (dashed), three (dot-dashed), and four (solid) loops in 
connection with Eq. (ESl) at the respective orders. 



Again it turns out to be useful to consider in addition to ( pTf ) the quantities where the 
scale invariant and the on-shell mass, respectively, has been used for the parameterization 
of the heavy quark: 
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The corresponding inverted relations read 
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As an example we compute m^^\Mz) for different number of loops and different match- 
ing points fith- The results can be found in Tab. |] where fith = Mz, Mf,, and 1 GeV has 
been chosen. It can clearly be seen that the four-loop result provides the most stable 
values for m^^\Mz)- 

A similar analysis as in Fig. |l] may be performed for the light-quark masses as well. 
For illustration, let us investigate how the /i'-^^ dependence of the relation between /ic = 
*^^'*(/ic) and m^^\Mz) changes under the inclusion of higher orders in evolution and 



m 



matching. As typical input parameters, we choose /ic = 1.2 GeV, Mb = 4.7 GeV, and 



a 



(Mz) = 0.118. We first evolve mf^{^) from /x = /ic to /i = /ij/i = ^^^> via Eq. (|T0|), 
then obtain m^^''(/i'^^)) from Eqs. (|3^) , and finally evolve m^^\jji) from /i = jj,^^^ to /i = Mz 

via Eq. (p^O]). In all steps, ai"^^(/i) is evaluated with the same values oi nj and as 

rri'^^\lj). In Fig. ^, we show the resulting values of mf^Mz) corresponding to A^-loop 
evolution with (A^ — l)-loop matching for A^ = 1, . . . , 4. Similarly to Fig. |l], we observe a 
rapid, monotonic convergence behaviour at the central matching scale /x^^-* = M^. Again, 
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number of loops 




(Mz) 




(running) 


{fith = Mz) 


ifith = Mb) 


ifith = 1 GeV) 


1 


0.6968 


0.7174 


0.7207 


2 


0.6025 


0.6051 


0.6236 


3 


0.5846 


0.5852 


0.5984 


4 


0.5798 


0.5801 


0.5783 



Table 5: Computation of m^\Mz) from jjc = fncil^c) = 1-2 GeV for different number of 
loops. For the matching scale Mz, Mb and 1 GeV has been chosen. 



the prediction for = 4 is remarkably stable under the variation of /x'-^^ as long as 
fx^^^ ^ 1 GeV. Fig. I can easily be reproduced with the help of the procedure mL2mH[] 
(see Appendix). 

If one chooses to perform the running of as{fi) with the help of A it is useful to have an 
equation at hand which relates this parameter in the full and effective theory. Combining 
Eqs. §), (I) and (|l]) one obtains § 

P',\n^ = {(3'^-(3,)l, + (b[-h)\nk-b[\n^ 

A Po 

+ Wib'i - h) In/, + 6f -bl-b', + b2 + l^ 
+ - &i) In' h + h,[-b[{b\ - 6i) 

+ &2 - ^2 - ^] In/. + i(-6f - - + 63) 

,,,,2 , ll^ 564731 82043^ 2633 1 

+ b\ (k + b'n-b2 ) H Cs ni } , (34) 

^ ^ 72^ 124416 27648^^ 31104 'j ^ ^ 

where /. = ln(/i^/A^) and the primed quantities refer to the {uf — l)-flavour effective 
theory. In this equation /x/j has been chosen for the matching scale which is particularly 
convenient, since it eliminates the renormalization group logarithms in (pi]). This choice 
is furthermore justified with the help of Figs. |1| and ^ where it can be seen that, in higher 
orders, the actual value of the matching scale does not matter as long as it is comparable 
to the heavy-quark mass. In Eq. (^) the four different powers in correspond to the 
different loop orders. Whereas at one-loop accuracy only the linear term in 1^ has to 
be taken into account at four- loop order also the 1//^ contribution has to be considered. 
Eq. ( p4D is implemented in the procedure DecLambdaDown [] . 

For completeness we also display the inverted relation of Eq. (^4|): 

/?oln^ = (/?o-« + (fei-&'i)ln/;,-6iln| 
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Figure 2: /i^^^ dependence of m^^\Mz) calculated from f^c = ^^c\l^c) = 1-2 GeV, Mb = 
4.7 GeV and af^Mz) using Eq. (||) at one (dotted), two (dashed), three (dot-dashed), 
and four (sohd) loops in connection with Eq. (^) at the respective orders. 
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(35) 



with = ln(/i|/(A')^). It is realized in the procedure DecLambdaUp [] . 

At this point we would like to mention that next to the coupling constant and quark 
masses also the gauge parameter and the quark and gluon fields obey decoupling relations. 
The corresponding equations and results can be found in [0. 
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5 Description of the main procedures 



In this section we describe the procedures which are most important for the practical 
apphcations, namely the combined running and decoupling of the strong coupling and 
the conversion of the on-shell mass to the MS one and vice versa. 

In RunDec.m some masses and couplings are set to default values which are used if 
they are not specified explicitly. They are collected in the set NumDef and read (also the 
corresponding symbol used in RunDec is given): 

Mtau : = 1.777 GcV , Mc : Mc = 1.6 GcV , Mb : Mi, = 4.7 GcV , 

Mt : Mt = 175 GcV , muc : /i^ = 1.2 GeV , mub : Hb = 3.97 GeV , (36) 

Mz: Mz = 91.18 GeV, asMz : ^(^^(Mz) = 0.118 . 

The following procedure computes a^J^\fi) where a^J^\fio) is used as input parameter. 
As input only agifJ^o), ^o, ^ and the number of loops have to be specified. Both n and 
m are determined according to the values of the quark masses given in NumDef. In case 
n ^ m the heavy quarks are consistently decoupled at the heavy quark scale itself where 
for the mass definition the on-shell scheme is used. 

• AsRunDec: 

— input: q;^"'^(//o), A*o, A*, number of loops 

— output: a^J^\ij) 

— uses: AlphasExact [] , A1L2A1H [] and A1H2A1L [] 

— comments: The decoupling is performed automatically at the pole mass of 
the heavy quark where the values defined in Numdef are taken. If /j, is lower 
than Mc, m = 3 is chosen, i.e. the strange quark is not decoupled. 

— example: In order to compute q;^^^(500 GeV) = 0.952 with four-loop accuracy 
if af\Mz) = 0.118 is given one has to use the command 

AsRunDec [asMz/ . NumDef , Mz/ . NumDef , 500 , 4] . 

The conversion of the on-shell mass M to the MS mass, m, can be computed with the 
help of the procedure mOS2mMS [] : 

• mOS2mMS: 

— input: M, Uf, number of loops 

— output: m^''f\M) 

— uses: AsRunDec [] andmOS2mMS[] (from the appendix) 

— comments: The relation is implemented up to order (three loops). Uf is 

the number of active flavours. as^^\M) is evaluated at the scale M where 
af'\Mz), as defined in Numdef, serves as a starting point. For the running 
and decoupling the procedure AsRunDec is used. 
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— example: In the case of the top quark the MS mass mt{Mt) = 164.6 GeV is 
obtained via mOS2mMS [175,6,3] where Mj = 175 GeV has been chosen. 

The inverted relation is implemented in 

• mMS2mOS: 

— input: Hm = fn^"'^Kl^m), nj, number of loops 

— output: M 

— uses: AsRunDec[] andmMS2mOS[] (from the appendix) 

— comments: The relation is implemented up to order (three loops), nf is 

the number of active flavours. a^r^\nm) is evaluated at the scale /i^ where 
af'^Mz), as defined in Numdef, serves as a starting point. For the running 
and decoupling the procedure AsRunDec is used. 

— example: In the case of the top quark the on-shell mass Mt = 174.7 GeV is 
obtained via mMS2mOS [165,6,3] where mt{Mt) = 165 GeV has been chosen. 

In the above procedures all quarks lighter than the one under consideration are as- 
sumed to be massless. More specialized procedures providing more freedom in the choice 
of parameters and the running presentations can be found in the Appendix. 
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Appendix: Detailed presentation of the Mathematica 
modules contained in RunDec 

In the following we list the procedures contained in the program package RunDec and 
provide a brief description. The order of the parameters specified in the field input cor- 
responds to the order required in the Mathematica procedures. The precision used for 
most of the numerical evaluations is controlled with the variable $NumPrec. For the proce- 
dures involving, e.g., numerical solutions of differential equations or recursive solutions of 
equations the default precision of Mathematica is kept which is for all practical purposes 
more than enough. Note that often the precision requested for with $NumPrec can not be 
reached when the input data are only known to a few digits. 

AsRunDec [] is not listed as it can already be found in Section ^. 
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Procedures related to the strong coupling constant 

• LamExpl: 

— input: as'^\lj)-, fJ^, rif, number of loops 

— output: A("/) 

— uses: Eq. 

— comments: — 

— example: From the knowledge of a^J^\Mz) = 0.118 the computation of 
A(^) = 0.2089 to three-loop accuracy proceeds as follows: 

LamExpl [asMz/ . NumDef , Mz/ . NumDef ,5,3]. 

• Lamlmpl: 

— input: a^J^^\fi), /i, Uf, number of loops 

— output: A("/) 

— uses: Eq. (||) 

— comments: Solves Eq. (|^) numerically for A'^"^^ 

— example: If af'\Mz) = 0.118 is given the computation of A*^^) = 0.2083 to 
three-loop accuracy proceeds as follows: 

Lamlmpl [asMz/ . NumDef , Mz/ . NumDef ,5,3]. 

• AlphasLam: 

— input: A^''^f\ fx, Uf, number of loops 

— output: af'^\ii) 

— uses: Eq. (H) 

— comments: An explicit warning is printed on the screen if the ratio fx/A^'^f^ is 
too small. 

— example: For A^^) = 0.208 and = 4.7 GeV the value of af^Mb) = 0.2163 
is obtained to three-loop accuracy with AlphasLam [0 . 208 ,4.7,5,3]. 

• AlphasExact: 

— input: ai"'^\fio), fiQ, /i, nj, number of loops 

— output: af'^\^) 

— uses: Eq. (|l|) 

— comments: Solves the differential equation numerically using Q;s(/io) as initial 
condition. An explicit warning is printed on the screen if the ratio /i/A^^f^ is 
too small where A^"-'') is obtained with the help of LamExpl [] . 
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— example: af'^M^) = 0.2167 is computed from a^f\Mz) through 
AlphasExact [asMz/ . NumDef , Mz/ . NumDef , Mb/ . NumDef ,5,3] where the 
three-loop formulae are used. 

Procedures relating different mass definitions 

• mOS2mMS: 

— input: M, {Mg}, a^J^^\n), fi, nj, number of loops 

— output: Tn^'^f\n) 

— uses: Eq. (|13[) and Tab. ^ 

— comments: The relation is implemented up to order (three loops). {Mg}. 
is a set of light quark masses which can also be empty. For consistency 
reasons their values must correspond to the on-shell mass. 

Note that the name of the procedure is the same as the one introduced in 
Section ^. The distinction is only in the number of the arguments. 

— example: The MS mass corresponding to the on-shell top quark mass of 
175 GeV is computed via mOS2mMS [175 , {} , . 107 , 175 , 6 , 3] where 
a^^-'(175 GeV) = 0.107 has been chosen. The result reads 

mt(175 GeV) = 164.64 GeV. Terms up to order have been used and light 
quark mass effects have been neglected. 

• mMS2mOS: 

— input: m^'^f\iJ,), {mg}, a^r^\n), fx, nj, number of loops 

— output: M 

— uses: Eq. (|17|) for general /i and Tab. ^ 

— comments: The relation is implemented up to order (three loops). In this 
case the light quark masses are defined in the MS scheme {mg}. 

Note that the name of the procedure is the same as the one introduced in 
Section ^. The distinction is only in the number of the arguments. 

— example: The on-shell mass corresponding to the MS top quark mass 

(175 GeV) = 165 GeV is computed via mMS2mOS [165 , {} , . 107 , 175 , 6 , 3] 
where a(6)(175 GeV) = 0.107 has been chosen. The result reads 
Mt = 175.35 GeV. Terms up to order have been used and light quark mass 
effects have been neglected. 

• mOS2mMSrun: 

— input: M, {Mg}, ak'^\ij), /i, Uf, number of loops 

— output: m^^f\n) 
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— uses: AlphasExact [] , mOS2mSI [] and mMS2mMS [] 

— comments: In a first step /i^ is computed and afterwards m^^f\n) is 
evaluated. The usage is identical to mOS2mMS [] . 

— example: (analog to mOS2mMS [] ) 
mMS2mOSrun: 

— input: m^^f\^), {mg}, as'^\^), /z, n/, number of loops 

— output: M 

— uses: AlphasLam [] , Lamlmpl [] , mMS2mMS [] and mMS2mOS [] 

— comments: In a first step /i^ is computed. Then Eq. (|T7D only has to be used 
for /i = /im- The usage is identical to mMS2mOS [] . 

— example: (analog to mMS2mOS [] ) 
mOS2mMSit: 

— input: M, {m^}, ai"^^(/i), n/, number of loops 

— output: m^'^^\iJi) 

— uses: Eq. (|T^ for general /i and Tab. ^ 

— comments: For the computation Eq. ([T7| ) is used in order to avoid the on-shell 
masses of the light quark masses {mq}. The usage is identical to mOS2mMS[]. 

— example: (analog to mOS2mMS []) . 
mOS2mSI: 

— input: M, {Mg}, ai"'^'*(M), n/, number of loops 

— output: fXrn = '^'^^^\l^m) 

— uses: Eq. (|I6|) and Tab. 

— comments: The scale invariant mass is computed from the on-shell mass. 

— example: In the case of the bottom quark, the mass = 3.97 GeV is 
evaluated via mOS2mSI [Mb/ . NumDef , {1 . 6} , . 217 , 5 , 3] where 
af'\Mi,) = 0.217 has been chosen. In the mass relation terms up to order 
have been used and quark mass effects arising from Mc/Mh with 

Mc = 1.6 GeV have been taken into account. 

mMS2mMS: 

— input: m^'^f\no), a^r^\iJ.o), a'r^\jj)^ Uf, number of loops 

— output: m^^f^fx) 

— uses: Eqs. (|) and 
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— comments: — 

— example: From mi,{Miy) = 3.85 GeV one finds mi){Mz) = 2.69 GeV with the 
help of mMS2mMS[3.85,0.217,asMz/.NumDef ,5,4] ) where af\M^) = 0.217 
and a^^\Mz) = 0.217 has been used. For the running the four-loop 
expressions have been used. 

• mMS2mSI: 

— input: m^'^f\fi), as'^\jj), jjL, nj, number of loops 

— output: Hm = ^ifJ'm) 

— uses: Eqs. (|) and (|10D 

— comments: The scale invariant mass is computed from the MS mass. 

— example: fib = 3.97 GeV is computed from the input uibiMb) = 3.85 GeV, 
Mb = 4.7 GeV and af'\Mb) = 0.217 via the command 

mMS2mSI [3 . 85 , . 217 , 4 . 7 , 5 , 4] . For the running the four-loop expressions 
have been used. 

• mMS2mRI: 

— input: m^^f\fi), nj, number of loops 

— output: m^^ 

— uses: inverted equation of ( pISD 

— comments: The relation is implemented up to order (three loops). 

— example: The regularization invariant mass, m^\Mz), corresponding to the 
MS bottom quark mass mfe(M^) = 2.695 GeV is computed via 
mMS2mRI[2.695,asMz/.NumDef ,5,3] where af)(Mz) =0.118 has been 
chosen. The result reads m^^Mz) = 2.872 GeV where terms up to order 
have been used. 

• mRI2mMS: 

— input: m^\jj)^ a^J^^\n), Uf, number of loops 

— output: m^^f \ij) 

— uses: Eq. (|18D 

— comments: The relation is implemented up to order (three loops). 

— example: The MS mass corresponding to the regularization invariant top 
quark mass of mf^ = 175 GeV is computed via 

mRI2mMS[175,0. 107, 175,6,3] where a ('^)( 175 GeV) = 0.107 has been 
chosen. The result reads mt(175 GeV) = 165.6 GeV where terms up to order 
ai have been used. 
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• mMS2mRGI: 



— input: m^'^f\jj,), a^' (/i), n/, number of loops 

— output: m 

— uses: Eq (0) 

— comments: — 

— example: The renormalization group invariant bottom quark mass 
corresponding to the MS mass mf'\Mz) = 2.69 GeV is computed via 
mMS2mRGI[2.69,asMz/.NumDef ,5,4] where af'^iMz) =0.118 has been 
chosen. The result reads m^*^^ = 14.25 GeV assuming four-loop accuracy. 

• mRGI2mMS: 

— input: m, Uf, number of loops 

— output: m^^f\^) 

— uses: Eq (|ll|) 

— comments: — 

— example: The MS mass corresponding to the renormalization group invariant 
bottom quark mass of 14.25 GeV is computed via 

mRGI2mMS[14.25,asMz/.NumDef ,5,4] where a'^^Mz) = 0.118 has been 
chosen. The result reads mi,{Mz) = 2.69 GeV assuming four-loop accuracy. 

Decoupling of the strong coupling and the masses 

At this point we once again want to stress, that the argument specifying the number 
of loops refers to the accompanied running, i.e. if "2" is chosen the decoupling relation 
is used to one-loop order. Furthermore, for the argument ruling the number of active 
flavours the number of light quarks, rii = nj — 1, is chosen. 

• DecAsUpOS: 

— input: af'\iith), Mth, fith, ni, number of loops 

— output: a'f^^^^fjith) 

— uses: Eq. (pSf ) 

— comments: For the heavy mass the on-shell definition is used. 

— example: The computation of af\Mz) = 0.1169 from the knowledge of 
af^Mz) =0.118 proceeds via 

DecAsUpOS[asMz/.NumDef ,175,Mz/.NuniDef ,5,4] where Mt = 175 GeV has 

been chosen and terms of order (indicated by the "4" in the last 
argument) have been included. 
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DecAsDownOS: 

- input: a'f-'^'^\fith) , Mth, fith, Ui, number of loops 

- output: afi\^th) 

- uses: Eq. (H) 

- comments: For the heavy mass the on-shell definition is used. 

- example: The computation of a^^''(200 GeV) = 0.1047 from the knowledge of 

) (200 GeV) = 0. 105 proceeds via DecAsDownOS [0 . 105 , 175 , 200 , 6 , 4] 
where Mt = 175 GeV has been chosen and terms of order (indicated by 
the "4" in the last argument) have been included. 

DecAsUpMS: 

- input: a'f-''^ (fith) , mhifJ'th), fJ'th, ni, number of loops 

- output: a'f^'~^'^\fith) 

- uses: Eq. (p3D 

- comments: The heavy mass is evaluated in the MS scheme at the scale ^th- 

- example: The computation of a^^\Mz) = 0.1170 from the knowledge of 
af\Mz) = 0.118 proceeds via 

DecAsUpMS[asMz/.NumDef ,165,Mz/.NumDef ,5,4] where mt{Mz) = 165 GeV 
has been chosen and terms of order (indicated by the "4" in the last 
argument) have been included. 

DecAsDownMS: 

- input: af''^'^\^th), mhi^J'th), fJ^th, ni, number of loops 

- output: a^s^'\fith) 

- uses: Eq. (ID 

- comments: The heavy mass is evaluated in the MS scheme at the scale ^th- 

- example: The computation of a^^''(200 GeV) = 0.1048 from the knowledge of 
af\200 GeV) = 0.105 proceeds via DecAsDownMS [0. 105, 165,200,6,4] 
where mt{Mz) = 165 GeV has been chosen and terms of order (indicated 
by the "4" in the last argument) have been included. 

DecAsUpSI: 

- input: a^'^'^fith), fJ^mth^ f^th, number of loops 

- output: af^+^\^ith) 

- uses: Eq. (HD 

- comments: Here the scale invariant mass /imt^ is chosen for heavy mass. 
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— example: (analog to DecAsUpMS [] ) 
DecAsDownSI: 

— input: a("'+^)(;Ut/i), /imt^, /ii/i, ni, number of loops 

— output: a^^'\iith) 

— uses: Eq. 

— comments: Here the scale invariant mass ^imth chosen for heavy mass. 

— example: (analog to DecAsDownMS [] ) 

DecMqUpOS: 

— input: m^^i\iith), a^f^'^fith), Mth, fJ^th, ni, number of loops 

— output: m^^^^^\jjith) 

— uses: Eq. ( p3|) 

— comments: For the heavy mass the on-shell definition is used. 

— example: The computation of m\ (Mz) = 2.697 GeV from 
mf'\Mz) = 2.7 GeV with order accuracy is performed via 
DecMqUpOS [2.7, asMz/ . NumDef ,175, Mz/ . NumDef ,5,4]. Here, 
af\Mz) = 0.118 and Mt = 175 GeV have been used. 

DecMqDownOS: 

— input: irS^'^^\^.th), oif''^^\^.th) , Mth, fJ'th, number of loops 

— output: m^^'-\fith) 

— uses: Eq. (H) 

— comments: For the heavy mass the on-shell definition is used. 

— example: The computation of ml^^Mz) = 0.583 GeV from 
m^^\Mz) = 0.58 GeV with order accuracy is performed via 
DecMqDownOS [0 . 58 , asMz/ . NumDef ,4.7, Mz/ . NumDef ,5,4]. Here, 
a(5)(M^) = 0.118 and Mb = 4.7 GeV have been used. 

DecMqUpMS: 

— input: m^^'\fith), oi^^'Kl^th), mhifJ'th), fJ'th, ni, number of loops 

— output: rn^^^^^\iith) 

— uses: Eq. ( PID 

— comments: The heavy mass is evaluated in the MS scheme at the scale ^ 

— example: (analog to DecMqUpOS [] ) 
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DecMqDownMS: 

— input: mf'^'^\iith) , af''^'^\iJith) , mhifJ'th), fJ'th, m, number of loops 

— output: mf^\iith) 

— uses: Eq. (p7|) 

— comments: The heavy mass is evaluated in the MS scheme at the scale ^th- 

— example: (analog to DecMqDownOS [] ) 

DecMqUpSI: 

— input: m^^'\iith), af'\iith), f^mth, f^th, Ui, number of loops 

— output: m''^'~^^\fith) 

— uses: Eq. (H) 

— comments: Here the scale invariant mass fimth chosen for heavy mass. 

— example: (analog to DecMqUpOS [] ) 

DecMqDownSI: 

— input: m^"'+^)(/it/,), a("'+^)(/it/,), ^mth, f^th, ni, number of loops 

— output: m'^^'\fith) 

— uses: Eq. (|9D 

— comments: Here the scale invariant mass fimth chosen for heavy mass. 

— example: (analog to DecMqUpMS [] ) 

DecLambdaUp: 

— input: A^"'''\ Hmth, "^u number of loops 

— output: A("'+i) 

— uses: Eq. (|§ 

— comments: For the heavy mass the scale invariant mass yUmt^ is used. 

— example: From A*^^) = 0.2876 one can compute A(^) = 0.208 with the help of 
DecLambdaUp [0.287, 3. 97, 4, 4] where = 3.97 and four-loop accuracy has 
been chosen. 

DecLambdaDown: 

— input: A*^"''^^^ IJ^mth-, ni, number of loops 

— output: A("') 

— uses: Eq. ( pSD 
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— comments: For the heavy mass the scale invariant mass /imt^ is used. 

— example: From A'-^-' = 0.208 one can compute A(^) = 0.208 with the help of 
DecLambdaDown[0.208,3.97,4,4] where fif, = 3.97 and four-loop accuracy 
has been chosen. 

Miscellaneous procedures 

The following modules provide some simple examples which mostly combine the modules 
described above. "L" stands for low and "H" for high. The condition I < h is assumed in 
all four procedures. 

• A1L2A1H: 

— input: a^sK^ii), /ii, {{n/^, Mt^^, /ij/^J, {n/j, Mt/^^, /itha}, • • •}, At2, number of 
loops 

— output: a^J^\iX2) 

— uses: AlphasExact [] and DecAsUpOS [] 

— comments: The set in the third argument may contain several triples 
indicating the number of flavours, the heavy (on-shell) quark mass and the 
scale at which the decoupling is performed. 

— examples: 1. For the computation of a^^''(500 GeV) = 0.0952 from 
a^'^^{Mc = 1.6 GeV) = 0.338 to 0{a^) accuracy the input would look as 
follows: A1L2A1H [0 . 338 , 1 . 6 , {{5 , 4 . 7 , 5} , {6 , 175 , 200}} , 500 , 4] Here, the 
matching is performed at 5 GeV and 200 GeV, respectively. 

2. Fig. |T| can be reproduced with the help of the following input 
AlL2AlH[0.36,1.777,{{5,4.7,mu5}},91.187,l] where / = 1,2,3,4 
corresponds to the number of loops and mu5/4.7 is the scale on the abscissa. 

• A1H2A1L: 

— input: af\^,i), fii, {{uf^, Mth^, fithA, {nf^, Mth^, fJ'thi}, ■ ■ ■}, /^2, number of 
loops 

— output: a^p(n2) 

— uses: AlphasExact [] and DecAsDownOS [] 

— comments: The set in the third argument may contain several triples 
indicating the number of flavours, the heavy (on-shell) quark mass and the 
scale at which the decoupling is performed. 

— example: Consider the inverse order of the first example of the previous 
procedure. The input 

A1H2A1L [0 . 0952 , 500 , {{6 , 175 , 200} , {5 , 4 . 7 , 5}} , 1 ■ 6 , 4] , indeed leads to 
a(^)(1.6 GeV) = 0.338. 
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• mL2mH: 

- input: mW(/ii), a^'^/ii), /ii, {{uf^, Mth,, fJ-thi}, {nf2, Mth^, /J-thi}, ■ ■ At2, 
number of loops 

- output: m^'*^(/i2) 

- uses: AlphasExact [] , mMS2mMS [] DecMqUpOS [] and DecAsUpOS [] 

- comments: The set in the fourth argument may contain several triples 
indicating the number of flavours, the heavy (on-shell) quark mass and the 
scale at which the decoupling is performed. 

- example: Using ai'^\l.2 GeV) = 0.403 and m(^)(1.2 GeV) = 1.2 GeV one 
finds m(5)(Mz) = 0.580 GeV with the help of 

mL2mH [1 . 2 , . 403 , 1 ■ 2 , {{5 , 4 . 7 , 5 . 0}} , Mz/ . NumDef , 4] . The decoupling of 
Mb = 4.7 GeV is performed at 5.0 GeV. In this way the results of Fig. ^ can 
be reproduced. 

• mH2mL: 

- input: mf\ni), af\ni), /xi, {{uf^, Mth^, fith^],{nf^, Mth^, fith^], . . /i2, 
number of loops 

- output: m^\ii2) 

- uses: AlphasExact [] , mMS2mMS [] DecMqDownOS [] and DecAsDownOS [] 

- comments: The set in the fourth argument may contain several triples 
indicating the number of flavours, the heavy (on-shell) quark mass and the 
scale at which the decoupling is performed. 

- example: Using q:(^)(M^) = 0.118 and m^^\Mz) = 0.580 GeV one finds 
m(^)(1.2 GeV) = 1.20 GeV with the help of 

mH2mL [0 . 580 , asMz/ . NumDef , Mz/ . NumDef , {{5 , 4 . 7 , 5 . 0}} , 1 . 2 , 4] . The 
decoupling of Mb = 4.7 GeV is performed at 5.0 GeV. 
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